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Abstract 



We propose a method to analyze infrared contributions to non-inclusive pro- 
cesses in QCD. We use the one-loop Sudakov form factor as a working exam- 
ple. Borrowing techniques from renormalization theory, we construct coun- 
terterms for the contributions from the soft and collinear regions, and we 
relate them to eikonal Wilson lines taken along non-lightlike directions. 
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I. INTRODUCTION 



Understanding the structure of hadronic final states in high energy colhsions involves 
applying QCD in situations with multiple mass scales. While the application of QCD to 
inclusive hard-scattering processes characterized by a single large mass scale is well under 
control, the use of QCD in less inclusive, multiple-scale processes is much subtler [|[|. 

A common feature of these more general situations is that soft gluon effects do not cancel. 
In tractable cases they merely factorize, and result in a more general factorization theorem, 
with new kinds of sensitivity to non-perturbative physics. An example is the transverse mo- 
mentum distribution for the Drell-Yan process 0. But we would like to treat, for example, 
Monte-Carlo event generators, where the theory does not yet go systematically beyond 
the leading logarithm approximation, despite the great practical importance of these event 
generators. The necessary extensions of the factorization theorems require that we devise 
new more precise methods. This is the aim of the present paper. 

We consider the well-known example of the one-loop part of the Sudakov form factor. 
Fig. [l|a. A full factorization theorem was obtained in [Q. Our method is a systematic 
subtraction method that generalizes the Bogoliubov /^-operation for renormalization. (For 
a related approach see Ref. [^].) The subtractions are obtained from matrix elements of 
certain gauge-invariant operators, which is essential if the factorization theorem is to be 
considered a real theorem of field theory. Moreover the subtractions are performed on the 
integrand, so that our final result for the ultraviolet part of the graph is infrared safe: it 
can be numerically integrated. Although our analysis is at the one loop level, the method 
should be generalizable to higher orders. 

The plan of the paper is as follows. In Sec. |T| we describe the method and the form of 
the counterterms, which are obtained from matrix elements of path-ordered exponentials of 
the gluon field along non-lightlike lines. In Sec. |T| we discuss the evolution due to changes 
in these directions. In Sec. 11^ we summarize and briefly discuss results and developments. 



II. THE INFRARED SUBTRACTION METHOD 

Consider the graph in Fig. |l|a. To simplify the calculations while retaining all the essential 
ingredients on which we want to focus, we work in a massive abelian theory with scalar 
quarks. We denote the quark mass by m and the gauge boson mass by m^. We work in a 
center-of-mass frame in which the incoming quark momenta pa and Pb are in the +z and 
—z directions, with 2p\p]^ = Q^. We will compute the asymptotic large Q behavior of the 
graph when the integrand is multiplied by a test function if{k/Q): 

P^^j ^ ^^2 f i2pA-k)-i2pB + k) ^^^^^^ 

■J (2 7r)4 (^k'^ _ ^2 _^ [(p^ _ /^)2 _ ^2 _^ ^ ^] [(^^^ _^ /^)2 _ ^2 _^ ^ 

— MS counterterm for UV divergence . (1) 

For the actual graph of Fig. |I]a we will have (f = 1. But retaining ip{k/Q) until the end 
of the calculation will assist us in identifying the nature of the terms we obtain. It will 
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FIG. 1. (a) Form factor graph; (b) soft and collinear regions. 



also enable us to immediately solve a more general class of problem, when an infrared finite 
vertex correction is inserted. 

Following the analysis of we identify the important regions of momentum space 

in r, as depicted in Fig. |T]b: 

1. soft region: all components of k'^ are much smaller than Q: X Q, where A is small; 

2. pA-collinear region: ~ Q, ~ A Q, k~ ~ A^ Q, with < k~^ < p\; 

3. ps-coUinear region: k~ Q, k± \ Q, k^ A^Q, with —p]^ < k^ <0; 

4. hard region: all components of k^ are of order Q: k^ ~ Q. 

We will construct a term for each region. We will do this starting from the smallest region 
(region 1, representing a 0-dimensional subspace) and systematically going to larger regions: 
regions 2 and 3 (1-dimensional) and region 4 (4-dimensional). The strategy for constructing 
the approximation for a given region will be 

a) take the original graph and subtract the terms for smaller regions; 

b) construct an approximation valid in the region being treated; 

c) show that it has no divergences or leading contributions from smaller regions; 

d) "renormalize" divergences associated with larger regions. 

Since, as we will see, only regions that are larger or smaller than the current region are 
relevant, the sum of all the terms will give a good leading-power approximation to the 
original graph. 
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A. Soft region 



We begin with the soft region. Expanding the integrand of T in powers of the small 
parameter A (in the definition of the region) gives 

2 11 

So = Zl^ f dk^dk-(fk— — — -ifi(0) . (2) 

The expansion in A gives the correct answer in the soft region unless one or both of the 
components k~^, k~ is too small. In that case k^k~ <^ k^, and only the quark propagators 
have significant dependence on k~^ and k~ . For the original graph, we deform the k^ and 
k~ contours out of the dangerous region. The approximation is only valid on the deformed 
contour, so that we must apply the same ze's to the approximation in Eq. (|^). Observe that 
the direction of the deformation of the contour integrations is determined by the fact that 
the interactions are initial-state interactions. Therefore the same deformation applies in all 
higher-order graphs as well. 

Now we observe that Eq. is the one-loop contribution to the vacuum expectation 
value of two Wilson lines taken along light-like directions pA, Pb- 

~ ^ {0\T[V,iPA)V,ipB)]\0) . . 

° {o\v,{PAm m^iPsm ' ^ ^ 

where self-energy contributions have been divided out and we have defined 

Vq{n) = V exp (^ig J dzA{zn)-n^ , Vq{n) = V exp (^—ig J dzA{zn)-n^ . (4) 

Our soft term 5*0 still has singularities from the collinear and ultraviolet regions. The 
ultraviolet singularity is to be dealt with by the standard subtractive renormalization pro- 
cedure, by adding a counterterm. We will now treat the collinear singularities in a similar 
fashion. The collinear singularities are entirely associated with the integration over the 
rapidity y = ln{k~^/k~) and are not removed by dimensional regularization [§]. 

As usual, the definition of the counterterms is non-unique. The reasons for our choice of 
the definition of the collinear counterterms are as follows. First, the new term is to have an 
operator interpretation, so that we can give it a nonperturbative definition. Secondly, the 
evolution equation is to be "simple" , in a sense to be specified later. Here we just observe 
that the new definition will involve considering Wilson lines along an arbitrary direction, 
and the evolution equation treats the dependence on this direction. 

Let us introduce two non-lightlike vectors ua, ub- ua = (uA^W^^Oj, ub = (u^,u]^,0^. 
Then we define the subtracted soft term by 

2 

S = I dk+ dk- rf^k ■ 



(2 7r)4 J (^k'^ - ml + i e) 



{k- -ie) {k+ + ie) 



Ub Ma 1 



^(0) 



{k~ — ie) {ub k+ + u%k~ + i e) {u\ k^ + k+ — ie) {k+ + ie) 
-UV . (5) 
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The first term in the square brackets is the unsubtracted soft term. The second term 
subtracts the divergence from the region colhnear to pa, i-e., from large positive rapidity y. 
The ub factor in this term provides a cut-off when the rapidity of k goes much below the 
rapidity of ub- Similarly, the third term subtracts the divergence from the region collinear 
to pb- Observe that, in order to avoid invalidating the contour deformation argument, we 
must have u\,Ub > 0, u'^,u'j^ < 0. Note that the collinear-to-pA subtraction term has no 
collinear-to-ps singularity; indeed it is power suppressed in this region. The same is true 
with A and B exchanged. Provided that we choose ua and ub to have large rapidities: 
''^a/'^^a ~ Pa/Pa^ ''^b/'^b ~ Pb/Pb^ then 5* continues to provide a good approximation 
to r in the soft region. The ultraviolet subtraction in the third line of Eq. (^) is understood 
to be done in the MS scheme. 

The cut-offs are gauge-invariant since, as we now show, they are obtained from a Wilson- 
line operators. Using the path-ordered exponentials introduced in Eq. (|^), we define 

^ ^ {o\v,iPA)vSBm (oiKMlo) (oi^Mlo) ... 

{0\V,{pA)Vg{uBm{0\V,{uA)V,{pB)\0) ' ^ ^ 

whose one-loop term gives Eq. (^). The first factor in the numerator gives the original soft 
term 5*0 in Eq. (^, and the two factors in the denominator, involving eikonal lines along 
spacelike directions, give the collinear subtractions. The remaining factors in the numerator 
cancel factors of a complete external propagator for the Wilson line. Given our one-loop 
formulae, this result appears to us to be unique, if we simply assume that the quantity which 
we calculate is the product of vacuum expectation values of some Wilson line operators. 



B. Collinear to A region 



Next consider the collinear-to-p^ region. We will construct a term A such that A + S 
gives a good approximation to F (valid to the leading power) over the whole region, including 
its soft subregion. It is obtained by expanding F — S in powers of the parameter A in the 
definition of the coUinear-to-p^ region, for which A;+ ~ Q, /c± ~ A Q, /c~ ~ A^ Q. Thus 



A 



(T-S) 
(2^ 



\pA-Co\l. 

dk+ dk~ S\i 



X 



(P - 



+ ie 



k+)2pB 



(-2p\k- + 2k+ [k- 



{k- - ie) {k+ + ie) 



m'^/{2p\)]-k'^ + ie) (2pBk+ + ie 

MO) ' 



{k~—ie) {uBk^ + u%k~ + ie) 



(p(fc+/Q,0,O) 
^(0) 



-uv 



(7) 



By this construction of A, we ensure that A-l- S* is a good approximation to the coUinear-to- 
Pa region. As to the soft region, observe that the first term in square brackets, the collinear 
expansion of the original graph, has a soft singularity at k~^ = 0. But the second term, 
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obtained from S, cancels this singularity. Notice that the third term is power suppressed in 
the soft region. Hence ^4+5* is a good approximation to the original graph in a neighborhood 
of both the soft and collinear-to-p^ regions. 

We now show that the A term, Eq. (|^), is power suppressed in region 3, where k is 
collinear to pb- This region merely overlaps with the region that this term is designed 
to treat. In the ps-collinear region, a good approximation is obtained by replacing the 
quark factor {2p\ — k~^)/[{pA_ — kY — rn? + ie)] by the light-like eikonal 2p^/ {—2p\k'^ — ie) = 
— 1/ {k~^+ie). This is the same replacement as for the soft region, so that we get a suppression. 
This is a result and solves what is effectively a problem in overlapping divergences: We 
obtained A by expanding about the line fc~ = k^- = — Fig. |I]b, and it is suppressed at 
the intersection, k = 0, with the collinear-to-ps line. The extrapolation back to the whole 
of the B region continues to be suppressed. 

We therefore have the following results: 

• The sum A + S gives a good approximation to F in regions 1 and 2. 

• The term A is power suppressed in regions that overlap with the A region or are smaller 
than it. 



C. Collinear to B region 

The expression for the ps-collinear term, B, is simply obtained from Eq. by inter- 
changing the indices -|- and — , A and B. 

Just as for the soft term S, the collinear terms can be obtained uniquely as the one-loop 
expansion of a product of gauge-invariant matrix elements: 

m,{PB) qmpA) {o\v,{PA)v,{uBm ^ _ mm v,iPA)\PB) (qikm v.ipBm 
m.ipA) VgiPBm {o\v,{uBm ' m.ipA) ysbM (oiv.mio) 

(8) 



D. Ultraviolet region 

Finally, we consider the hard part of the form factor. We first observe that the sum of 
the terms for the smaller regions, A + B + S, provides a good approximation to F in the 
whole of the smaller regions 1, 2, and 3. To obtain a term H for the hard region (4), we 
simply take the massless approximation toT — A — B — S 
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Note that the all non-lightlike eikonal terms have canceled. 

To verify that the subtractions are correct, we perform the k^, k 
case This gives 

1 -kVQ^ \^ ri + R 



— V^(A;+/Q,0,O) 
ie) 



(9) 



integrations in the 
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if 4kVg2 > 1 



(10) 



(11) 



/4kVg2 _ ; 

Observe that, as a result of the infrared subtractions, the k^ integration in Eq. ([To|) is regular 
at small k^. The large k^ behavior, on the other hand, has to be dealt with via an ultraviolet 
subtraction. 

Eqs. (|^), (|^ and provide the main result of the application of our method. The 
analysis of Refs. 0,^,0] establishes the list of leading regions. For each region, we have 
defined a term that, added to the terms for smaller regions, gives a good leading-power 
approximation to the original graph in that region. For each term, we have shown that 
it does not receive leading contributions from regions that are smaller or have an overlap 
with the region being considered. Thus the sum of all the terms gives a good leading-power 
approximation to the original graph: 

H + A + B + S = T + nonleading power . (12) 



III. EVOLUTION 

The spacelike vectors ua, ub in the soft and coUinear subtractions are not physical 
parameters, and the sum H + A + B + S is independent of them. The utility of these vectors 
comes from the fact that evolution equations in i/a = \\ti\u\/u'X\ and i/b = |ln|M^/M^| 
can be applied to the factors in the factorization formula to extract effects associated with 
large logarithms [^Q. For example, differentiating Eq. (H) with respect to i/a gives 

dS ^ ig^ r d^k 

dvA {2tiY J {k^ -ml + ie) {u\k- + UAk+ -ief 

This is to be contrasted with the evolution equations obtained in , which suffer from power- 
law corrections. Our new method has all the power law corrections in the initial construction 
of the A, B, S and H terms. Since no further approximations are needed in solving the 
evolution equations 0, this can make for easier and cleaner proofs in more complicated 
situations. It also makes it easier to locate the sources of the power-law corrections: they 
are all associated with the approximations made to derive the factorization. 
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IV. OUTLOOK 



In this paper we have proposed a systematic subtraction procedure to analyze infrared 
contributions to non- inclusive scattering processes in QCD. The counterterms arise from 
gauge-invariant Wilson line operators which are in spacelike directions. The form of these 
operators allows simple evolution equations to be obtained, which have no power-law cor- 
rections; this makes it much easier to pinpoint where the power-law errors in factorization 
theorems occur. Observe that, even in our one-loop example, the pattern of infrared contri- 
butions in the coUinear regions is equivalent to that of overlapping divergences in renormal- 
ization theory. 

Uniqueness of the operator formulae follows from the assumption of a factorized form 
for the complete amplitudes. Given the general and systematic nature of our approach, 
we expect to generalize it to all higher orders and to other, more physically interesting, 
processes. 

This research is supported in part by the US Department of Energy under grant No. DE- 
FG02-90ER-40577. 
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